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-y^ . Abstract 

PLh ■ A coupling by reflection of a time-inhomogeneous diffusion process on a manifold 

are studied. The condition we assume is a natural time-inhomogeneous extension of 
lower Ricci curvature bounds. In particular, it includes the case of backward Ricci 
flow. As in time-homogeneous cases, our coupling provides a gradient estimate of 
the diffusion semigroup which yields the strong Feller property. To construct the 
coupling via discrete approximation, we establish the convergence in law of geodesic 
J> ■ random walks as well as a uniform non-explosion type estimate. 

in 

CN ■ 1 Introduction 

q 

In stochastic analysis, coupling methods of stochastic processes have played a prominent 
role in the literature. Given two stochastic processes Yx(t) and Y 2 (t) on a state space M, 
a coupling X(£) = (Xi(t), X 2 (t)) of Yi(t) and l^(t) is a stochastic process onMxM such 
that Xi has the same law as Yi for i — 1,2. By constructing a suitable coupling which 
reflects the geometry of the underlying structure, one can obtain various estimates for 
heat kernels, harmonic maps, eigenvalues etc. under natural geometric assumptions (see 
[11, 14, 27] for instance). Recently, the heat equation on time-inhomogeneous spaces such 
as Ricci flow have been studied intensively (see [1, 7, 18, 17, 19, 24, 29] and references 
therein). These studies have succeeded in revealing a tighter connection between the heat 
equation and the underlying geometric structure even in time-inhomogeneous cases. It 
should be remarked that an idea of coupling methods lies behind some of them [1, 17, 19, 
24] in an essential way . 

This paper is aimed at constructing a coupling by reflection of a diffusion process 
associated with a time-dependent family of metrics such as (backward) Ricci flow. Let 
M be a smooth manifold with a family of complete Riemannian metrics {<7(£)}tepi,T 2 ]- 
By {-X"(£)}tem,T 2 ]> we denote the g(£)-Brownian motion. It means that X(t) is the time- 
inhomogeneous diffusion process on M associated with A 9 ( t )/2, where A s ( t ) is the Lapla- 
cian with respect to g(t) (see [7] for a construction of g(£)-Brownian motion). As in 
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time-homogeneous cases studied in [8, 13, 16, 25, 26, 27] under a lower Ricci curvature 
bound, a coupling by reflection X(t) = (X 1 (t),X 2 (t)) of two g(£)-Brownian motions start- 
ing from a different point provides us a useful control of the coupling time r*, the first 
time when coupled particles meet. A simple version of our main theorem which states 
such a control is as follows: 

Theorem 1.1 Suppose 

d t g(t) < Ric 9(t) (1.1) 

holds. Then, for each x±,x 2 G M, there exists a coupling X(t) := (Xi(t),X 2 (t)) of two 
g(t)-Brownian motions starting at (xi,x 2 ) satisfying 



P[r* > t] < P 



' inf B{) _<W*i,* 2 ) 



1.2) 



for each t, where dg^) is the distance function on M with respect to g{T\) and B(t) is a 
1 -dimensional standard Brownian motion starting at the time T\. 

For the complete statement of our main theorem, see Theorem 4.1. There we also study 
a diffusion process which generalizes the g(£)-Brownian motion. The condition (1.1) can 
be interpreted as a time-inhomogeneous analogue of nonnegative Ricci curvature (see 
Remark 4.2). This condition is essentially the same as backward super Ricci flow in [19] 
(Our condition is slightly different in constant since our g(£)-Brownian motion and hence 
the heat equation corresponds to A 9 ( t )/2 instead of A 9 ( t )). Obviously, (1.1) is satisfied if 
g(t) evolves according to backward Ricci flow dtg{t) = Ric 9 (t). As in the time-homogeneous 
case, our coupling time estimate yields a gradient estimate of the heat semigroup which 
implies the strong Feller property for the heat semigroup (see Corollary 4.3). Note that, 
when g(t) is a backward Ricci flow, the same estimate as Corollary 4.3 is also obtained in 
[7] by using techniques in stochastic differential geometry. 

To explain our approach to Theorem 1.1, let us review a heuristic idea of the construc- 
tion of a coupling by reflection as well as that of the derivation of (1.2). Given a Brownian 
particle X\, we will construct X 2 by determining its infinitesimal motion dX 2 (t) G Tx 2 (t)M 
by using dXi(t) G T Xl (t)M. First we take a minimal g(£)-geodesic 7 joining Xi(t) and 
X 2 (t). Next, by using the parallel transport along 7 associated with the g(£)-Levi-Civita 
connection, we bring dX\{t) into Tx 2 {t)M. Finally we define dX 2 (t) as a reflection of it 
with respect to a hyperplane being g(t)-perpendicular to 7 in Tx 2 (t)M. From this con- 
struction, the Ito formula implies that d g u){Xi(t), X 2 (t)) should become a semimartingale 
at least until (X 1 (t),X 2 (t)) hits the g(£)-cutlocus Cut 9 ( t ). The semimartingale decompo- 
sition is given by variational formulas of arclength. On the bounded variation part, there 
appear the time-derivative of d g <t) and the second variation of d g m , which is dominated in 
terms of Ricci curvature. With the aid of our condition (1.1), these two terms are com- 
pensated and a nice domination of the bounded variation part follows. Thus the hitting 
time to of d g (t)(Xi(t), X 2 (t)), which is the same as r*, can be estimated by that of the 
dominating semimartingale. Indeed, we can regard 2B(t) + (i 9 (Ti)(^i,a ; 2) which appeared 
in the right hand side of (1.2) as the dominating semimartingale. The effect of our re- 
flection appears in the martingale part 2B(t) which makes it possible for the dominating 
martingale to hit 0. This construction seems to work as long as (Xi{t),X 2 (t)) is not in 
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the cutlocus. Moreover, we can hope it possible to construct it beyond the cutlocus since 
the cutlocus is sufficiently small and the effect of singularity at the cutlocus should make 
d g ^(Xi(t) , X 2 (t)) to decrease. If we succeed in doing so, the bounded variation part will 
involve a "local time at Cut g ( t )". It will be negligible since it would be nonpositive. We 
can conclude that almost all technical difficulties are concentrated on the treatment of 
singularity at the cutlocus in order to make this heuristic argument rigorous. In fact, 
Theorem 1.1 is shown in [21] by using SDE methods when the g(£)-cutlocus is empty for 
every t G [Ti,T 2 ]. 

Our construction of a coupling by reflection is based on a time-discretized approxima- 
tion as studied in [16, 25]. We construct a coupling of geodesic random walks each of whose 
marginals approximates the original diffusion process. The construction will be finished 
after taking a limit so that these approximations converge. Our method has a remarkable 
advantage in treating singularities arising from the cutlocus. In our construction, we can 
avoid to extract a local time at Cut s (*) and directly obtain a dominating process which 
does not involve such a term. In the present framework, the singular set Cut fl ( t ) also de- 
pends on time parameter t and hence treating it by using stochastic differential equations 
seems to be more complicated than in the time-homogeneous case. 

Different kinds of couplings are studied in above-mentioned papers. Based on the 
theory of optimal transportation, McCann and Topping [19, 24] studied a coupling of heat 
distributions which minimizes their transportation cost. They used the squared distance in 
[19] or Perelman's C- functional in [24] respectively to quantify a transportation cost. Their 
coupling is closely related to coupling of Brownian motions by parallel transport along 
minimal (£-)geodesics. In fact, studying a coupling by parallel transport by probabilistic 
methods recovered and extended (a part of) their results in [1] and [17] respectively. Note 
that our approach via time-discretized approximation is used in [17]. In addition, we also 
can construct a coupling by parallel transport by using our method to recover a result in 
[1] (see Theorem 4.6). It explains that our approach is also effective even when we study 
a different kind of couplings. 

We give a remark on a difference in methods between ours and Arnaudon, Coulibaly 
and Thalmaier's one [1] to construct a coupling by parallel transport. They consider one- 
parameter family of coupled particles along a curve. Intuitively saying, they concatenate 
coupled particles along a curve by iteration of making a coupling by parallel transport. 
Since "adjacent" particles are infinitesimally close to each other, we can ignore singular- 
ities on the cutlocus when we construct a coupled particle from an "adjacent" one. It 
should be noted that their method does not seem to be able to be applied directly in order 
to construct a coupling by reflection. Indeed, their construction of a chain of coupled par- 
ticles heavily relies on a multiplicative (or semigroup) property of the parallel transport. 
However, our reflection operation obviously fails to possess such a multiplicative property. 
Since our reflection map changes orientation, there is no chance to interpolate it with a 
continuous family of isometries. 

In what follows, we will state the organization of this paper. In the next section, we 
show basic properties of a family of Riemannian manifolds ((M,g(t)))t- In particular, 
we prove that Riemannian metrics {g(t)) t are locally comparable with each other. It will 
be used to give a uniform control of several error terms which appear as a result of our 
discrete approximation. In section 3, we will study geodesic random walks in our time- 
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inhomogeneous framework. There we introduce them and prove the convergence in law to 
a diffusion process. After a small discussion at the beginning of the section, the proof is 
divided into two main parts. In the first part, we will give a uniform estimate for the exit 
time from a big compact set of geodesic random walks. Our assumption here is almost the 
same as in [18] where non-explosion of the diffusion process is studied (see Remark 3.3 (ii) 
for more details). In the second part, we prove tightness of geodesic random walks on the 
basis of the result in the first part. In section 4, we will construct a coupling by reflection 
and show an estimate of coupling time, which completes the proof of Theorem 1.1 as a 
special case. 

2 Properties on time-dependent metric 

Let M be a m-dimensional manifold. For — oo < T\ < T 2 < oo, Let (<?(0)te[Ti,T 2 ] be a 
family of complete Riemannian metrics on M which smoothly depends on t. 

Remark 2.1 It seems to be restrictive that our time parameter only runs over the compact 
interval [T , 1 ,T 2 ]. An example of g(t) we have in mind is a solution to the backward Ricci 
flow equation. In this case, we can work on a semi-infinite interval [T l5 oo) only when 
we study an ancient solution of the Ricci flow. Thus T 2 < oo is not so restrictive. In 
addition, we could extend our results to the case on [Ti,oo) with a small modification of 
our arguments. It would be helpful to study an ancient solution. To deal with a singularity 
of Ricci flow, it could be nice to work on a semi-open interval (T , 1 ,T 2 ] ; where T\ is the 
first time when a singularity emerges. In that case, we should be more careful since we 
cannot give "an initial condition at T\" to define a g(t)-Brownian motion on M. 



We collect some notations which will be used in the sequel. Throughout this paper, 
we fix a reference point o G M. Let No be nonnegative integers. For a, b G 1R, a A b 
and a V b stand for min{a, b} and max{a, b} respectively. Let Cut 9 ( t )(:r) be the set of the 
g(£)-cutlocus of x on M. Similarly, the g(£)-cutlocus Cut fl ( t ) and the space-time cutlocus 
Cut st are defined by 

Cut fl(t ) := {(x,y) G M x M \ y G Cut ff(t) (x)} , 
Cutsx := {(t,x,y) G [T U T 2 ] x M x M \ (x,y) G Cut fl(t) } . 

Set D(M) := {(x, x) \ x G M}. The distance function with respect to g(t) is denoted by 
d g (t)(x,y)- Note that CutsT is closed and that g? 9 (-)(-, •) is smooth on [Ti,T 2 ] x M x M \ 
(Cut ST U[Ti,T 2 ] x D{M)) (see [19], cf. [18]). We denote an open #(s)-ball of radius R 
centered at x G M by B^ (x). Some additional notations will be given at the beginning 
of the next section. 

In the following three lemmas (Lemma 2.2-Lemma 2.4), we discuss a local comparison 
between d g (t) and d g ^ for s^t Those will be a geometric basis of the further arguments. 

Lemma 2.2 Let M be a compact subset of M. Then there exists k, = k(M ) such that 
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holds on M for t,s G [Ti , T 2 ] . In particular, if a minimal g(s) -geodesic 7 joining x, y G M 
zs included in M 0; £/ien ; /or £ G [Ti,T 2 ], 

d g (t)(x,y) < e Klt - sl d g{s) (x,y). 
Proof. Let 7r : TM — > M be a canonical projection. Let us define M by 
M := {(t,v) G [T U T 2 ] x TM | tt(u) G M , \v\ g{t) < l}. 

Note that M is closed since (?(•) is continuous. We claim that M is sequentially compact. 
Let us take a sequence ((t n ,v n )) n &n C M . We may assume t n — > t G [Ti,T 2 ] and 
7r(f n) ->?£ Mo as n — > 00 by taking a subsequence if necessary. Let U be a neighborhood 
of p such that {r> G TM | 7r(i>) G U} ~ [/ x IR"\ For sufficiently large n, we regard r n 
as an element of U X IR m and write r n = (p n ,v n ). If we cannot take any convergent 
subsequence of (v n ) n£ ^, then \v n \ — >■ 00 as n — > 00, where | • | stands for the standard 
Euclidean norm on M. m (irrelevant to (<7(£))te[7i,r 2 ])- Set t^ = (p n , \v n \~ l v n ). Then, there 
exists a subsequence (i^JfceN C (f'J ne N such that i>' — >■ f^, = (p, i>') as n — > 00 for some 
v' G M m with |r'| = 1. Since g(-) is continuous, g(t nk )(v' nk ,v' ) — > gifyiv'^^v'^) as k — >■ 00. 
On the other hand, g{t nk ){ v ' nk -> v 'n k ) < |£>n fc | -2 — > since g(t n )(v n ,v n ) < 1. Thus r 7 must 
be 0. It contradicts with \v'\ = 1. Hence Mo is sequentially compact. 

Since M E5 (t,v) h-» d t g(t)(v,v) is continuous, there exists a constant k = k(M ) > 
such that |<9 t (?(£)(i>, v)\ < 2k for every (t,v) G M . Take r G 7r -1 (M ), i> 7^ 0^). Then 

d t g(t)(v,v) = \v\ 2 g{t) d t g(t){\v\-( t) v,\v\-( t) v) < 2«|v|* (t) . 

Thus dtlogg(t)(v,v) < 2k holds. By integrating it from s to i with s < t, we obtain 
g(t)(v,v) < e 2Kiyt ~ s ^g(s)(v ,v ). We can obtain the other inequality similarly. 
For the latter assertion, for a, b with 7(a) = x and 7(6) = y, 

rb rb 

dg(t)(x,y)< / \i(u)\ (t) du<e Klt ~ sl \^(u)\ q(s) du = e Klt -^d g(s) (x,y). 



J a J a 

D 

Lemma 2.3 For i? > ; x G M and £ G [Ti, T 2 ] ; £/iere exists 5 = 5(x, t, R) > snc/i i/W 
T^ s) (x) C £3? Or) forr < R and s G [Ti,r 2 ] wjf/i |s - t\ < 5. 

Proof. Set k := k(B^(x)) as in Lemma 2.2 and 0" := /t _1 log2. Take p G Br (x) 
and a minimal g(s)-geodesic 7 : [a, b] — > M joining x and p. Suppose that there exists 
Mo G [a,b] such that 7(1x0) £ -^L ( x ) c - Let w := inf{w G [a, 6] | 7(w) G -B 3r (x) c }. Since 
7([a, «o]) C B^J(x) C -83^(2;) and d g (t)(x, 7(1*0)) — 3r, Lemma 2.2 yields 

wo ruo q„ 



d g{s) (x,p)> \j(u)\ g(s) du> e K / |7(w)| ff(t) dw=y. 

This is absurd. Hence 7([o, b]) G B^J(x). In particular, 7(6) = p G S 3r .(x). D 
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Lemma 2.4 For R > 0, there exists a compact subset M = M (R) of M such that 



pe M 



te 



mi d g{t) (o,p)<R\cM . 



(2.1) 



Proof. For each t G [T^TVIj take S(o,t,R+ 1) > according to Lemma 2.3. Take 
{*i}jLi C [Ti,T 2 ] such that 

n 

[T 1} T 2 ] c \J(ti - 5{o, U, R + 1),U + 5{o, U, R + l)). 
i=i 

Let us define a compact set M C M by M := UILi -^3-R (°)- Take p £ M such that 
inf Tl < t < Ta d g (t)(o,p) < R. For e G (0,1), take s G [7i,T 2 ] such that d g ^(o,p) < R + e. 
Then there exists j G {1, . . . , N} such that \s — tj\ < S(o, tj, R + 1). By Lemma 2.3, it 

implies p G B^ £ (o) C B^ +e Jo) C UI=i ^3(r+ £ )( )- Hence the conclusion follows by 
letting £ 4 0. □ 

Another useful consequence of Lemma 2.2 and Lemma 2.3 is the following: 
Lemma 2.5 d g (.\(-, •) is continuous on [Ti,T2] x M x M. 

Proof. Since the topology on [Ti,T 2 ] x M x M is metrizable, It suffices to show 
lim^oo d g( t n )(x n ,y n ) = d g{t) (x,y) when (t n ,x n ,y n ) ->■ (i,x,y) as n ->• oo. By the tri- 
angle inequality, 

Take i? > so that f?y(a;) includes a minimal g(£)-geodesic joining x and y. Take k = 
n^B^iyx)) according to Lemma 2.2. We can easily see that every minimal g(£)-geodesic 
joining y and y n is included in B^^x) for sufficiently large n G N. Thus Lemma 2.2 yields 

limsupd ff(tn )(y,y n ) < limsupe K|t ~ tn| d 9 ( t) (?/,?/ n ) = 0. 

n—too n— ►oo 

We can show d g u n )(x, x n ) —¥ similarly. Take a minimal g(£ n )-geodesic 7„ : [a, b] — > M 
joining x and y. By our choice of R, Lemma 2.2 again yields 

d 9 (t n )(x, ln {u)) < d g{tn) (x,y) < e K l*-*"lrf 9(t) (x,y) < e^'^R. 

It implies limsup n ^ 00 (i 9 (t n )(a;, y) < d g ( t )(x,y). In addition, 7„ is included in B\^L(x) for 

sufficiently large n. Thus Lemma 2.3 and Lemma 2.2 yield d g ( t )(x,y) < e K ^ tn ^d g ^ tn )(x,y). 
Hence the conclusion follows by combining these estimates with (2.2). □ 

Before closing this section, we will provide a local lower bound of injectivity radius which 
is uniform in time parameter. 

Lemma 2.6 For every Mi C M compact, there is f = f (Mi) > such that d g ( t )(y, z ) < 
f implies (t,y,z) ^ CutsT f or any (t,y,z) G [Ti,T2J x Mi x Mi. 
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Proof. Take R > 1 so that sup iG [ TlT2 ] sup xgMl d g (t)(o,x) < R — 1. By Lemma 2.4, there 
exists a compact set M C M such that (2.1) holds. For every t G [7\,T 2 ] and x G M 1; 
(t,x,x) <£ CutsT- It implies that there is r) tx G (0, 1) such that (s,y,z) <£ CutsT whenever 
d g (t)(x,y) V d g (t)(x,z) V \t — s\ < r) tjX since Cutsx is closed. Thus there exist N G N and 
(ti,Xi) G [T!,^] X Mi (z = 1, . . . , JV) such that 



TV 
i=l 



[Ti,T 2 ] X Mi C |J (ii - ^p,t, + ^p) x b£" ,,(.r, 



Set r > by 

1 / « 



r n := - exp max m. T . mm m. _., 

2 V 2 i<*<Af " */ i<i<v *' * 



where k = «(Mo) > is as in Lemma 2.2. Take (s,y,z) G [Ti,T 2 ] x Mi x Mi with 
dg(s)(y,z) < f . Take j G {1,...,N} so that \s - tj\ V d g ( ti )(xj,y) < %,^/2. By virtue 
of the choice of i? and M , Lemma 2.4 yields that every g(s)-geodesic joining y and z is 
included in M . Thus Lemma 2.2 yields 



d g ( tj )(y,z) <e Kls ^d g{s) (y,z) < 



r lt, 



2 
It implies |s - tj\ V d g ( tj )(xj,y) V d g ( tj )(xj, z) < 7] tj , Xj and hence (s, y, 2) ^ Cut S T- □ 

3 Approximation via geodesic random walks 

Let (Z(t))t£[T lt T 2 ] b e a family of smooth vector fields continuously depending on the pa- 
rameter t G [Ti, T 2 ]. Let X(t) be the diffusion process associated with the time-dependent 
generator Jf t = A g r t -\/2 + Z(t) (see [7] for a construction of X(t) by solving a SDE on the 
frame bundle). Note that (t,X(t)) is a unique solution to the martingale problem associ- 
ated with dt + Jif. on [Ti,T 2 ] x M (see [11] for the time-homogeneous case. Its extension 
to time-inhomogeneous case is straightforward; see [23] also). 

In what follows, we will use several notions in Riemannian geometry such as exponen- 
tial map exp, Levi-Civita connection V, Ricci curvature Ric etc. To clarify the dependency 
on the metric g(t), we put (t) on superscript or g(t) on subscript. For instance, we use 
the following symbols: exp™, V™ and Ric 9 (t). We refer to [6] for basics in Riemannian 
geometry which will be used in this paper. 

For each t G [T\,T 2 ], we fix a measurable section $ (t) : M -» <^ (t) (M) of the g(t)- 
orthonormal frame bundle ff^\M) of M. Take a sequence of independent, identically 
distributed random variables {£„}„ e N which are uniformly distributed on the unit disk 
in R m . Given xq G M, let us define a continuously-interpolated geodesic random walk 
(X a (t))te[T 1 ,T 2 ] on M starting from xq with a scale parameter a > inductively. Let 
tt ] := (Ti + a 2 n) A T 2 for n G N . For t = T x = t { *\ set X a {T 1 ) := x . after X a (t) is 
defined for t G [T 1: tn ], we extend it to t G [£« , £„+i] by 

|„ +1 := v / ST24<' i "''(.Y'>(«if ) ))?„ + i, 
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For later use, we define N^ := inf {n G No | t n+ i — t n < a 2 }. This is the total number of 
discrete steps of our geodesic random walks with scale parameter a. Set ff := C([Ti, T 2 ] — > 
M) and <2) := D([Ti,T 2 ] —¥ M). By using a distance d 9 m) on M, we metrize ^ and Ql 
as usual so that ^ and Q) become Polish spaces (see [9] for a distance function on 3>, for 
example). Set ^ := C([T 1 ,T 2 ] — > [0, oo)). Let us define a time-dependent (0, 2)-tensor 
field {vZ(t)f by 

(VZ(t))\X,Y) := \ ((V$Z(t),X) g{t) + (W?Z(t),Y) g{t) 



Assumption 1 There exists a locally bounded nonnegative measurable function b on 
[0, oo) such that 

(i) For all t G [T^^), 

2(VZ(t)) b + d t g(t) < Ric 9(t) +b(d g{t) (o, .))g(t). 

(ii) For each C > 0, a 1- dimensional diffusion process y t given by 

dy t = d(3 t + - (c + / b(s)ds ) dt 

does not explode. (This is the case if and only if 

r°° / rv \ f y ( f z \ 

exp I — / h(z)dz J / exp I / h(£)dt; J dzdy = oo, 

where b(y) := C + j^b(s)ds. see e.g. [12, Theorem VI. 3.2].) 
Our goal in this section is to prove the following: 

Theorem 3.1 Under Assumption 1, X a converges in law to X in ^ as a — > 0. 
Most of arguments in this section will be devoted to show the tightness i.e. 
Proposition 3.2 (X a ) Q , e ( 0i i) is tight in "H? . 
In fact, as we will see in the following, Proposition 3.2 easily implies Theorem 3.1. 

Proof of Theorem 3.1. By virtue of Proposition 3.2, for any subsequence of (X a ) Q , e ( 0) i) 

there exists a further subsequence (X afc )fc e jj which converges in law in c io as k — > oo. Thus 
it suffices to show that this limit has the same law as X. Let (/3 Q (t)) ie [o,oo) be a Poisson 
process of intensity a~ 2 which is independent of {£ n }neN- Set 

P a (t):=(T 1 + a 2 (3 a (t-T 1 ))At^i ) . 
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Then the Poisson subordination X ak ((3 ak (-)) also converges in law in Qi to the same limit 
(see [4] for instance). Note that (j3 a (t),X a (J3 a (t))) te < TliT2 ] is a space-time Markov process. 
The associated semigroup P t and its generator J2? ^ are given by 

^) /:=e - ( ^^((LM (g (.)y /+ £ ((t "y a " 2) ' ( g (a y (a 7y 

& a) f:=a- 2 (qMf-f), 

where 

g ( Q )/(t, x) := E [/(t + a 2 , expj? (a«s/ro + 2$W (x)£i + a 2 Z(t)) . 

We can easily prove JgfW/ -»■ (fy + Sf.)f uniformly as a -»■ for / G Cg°([7i,T2] x 
M). Since (/3 Q (t),X a ( / 9 a (t))) tG [T 1 ,T 2 ] is a solution to the martingale problem associated 
with Jzf^, the limit in law of (/3 ak (t),X a (/3 ak (t)))t£[T 1 ,T 2 } solves the martingale problem 
associated with <9 t + Jz^. By the uniqueness of the martingale problem, this limit has the 
same law as that of (t,X(t))t£[T lt T2]- It completes the proof. □ 

Remark 3.3 (i) A result on a convergence of semigroups [15] was used to show the 
convergence of finite dimensional distributions in the time-homogeneous case [5] (see 
[25] also). It is not so clear that we can employ the same argument in our time- 
inhomogeneous case. One difficulty arises from the absence of invariant measures 
for semigroups even in the case Z(t) = 0. Although the g(t)-Riemannian measure is 
a unique invariant measure for A g ^, this measure also depends on time parameter. 
Thus we cannot expect that it becomes an invariant measure of semigroups. This 
obstacle also prevents us to employ the existing theory of time- dependent Dirichlet 
forms (see [20] for instance) in order to study our problem. 

(ii) Proposition 3.2 also asserts that any subsequential limit in law is a probability mea- 
sure on *€ . Since we have not added any cemetery point to M in the definition of 
'io ' , Theorem 3.1 implies that X cannot explode. It almost recovers the result in [18]. 
Our assumption is slightly stronger than that in [18] on the point where we require 
(ii) for all C > 0, not a given constant. Note that we will use Assumption 1 (ii) 
only for a specified constant 2Co given in Lemma 3. 9. However, its expression looks 
complicated and it seems to be less interesting to provide a explicit bound. 

Now we introduce some additional notations which will be used in the rest of this paper. 
For t G [Ti,T 2 ], we define [t\ a by \t\ a := sup {a 2 n + T\ \ n G No,« 2 n + 7\ < t}. Set 
& n ■= o-(fi, . . . i n ). For R > 1, let us define a R : % ->■ [T U T 2 ] U {00} by 

a R (w) := inf {t G [T U T 2 ] \ w(t) > R - 1} , 

where inf = 00. We write Or := a R {d g ^{o, X a {-))) and a R := a^ 2 ([a R \ a — T x ) + 1. 
Note that a R is an #"„-stopping time. For each t G [Ti,T 2 ] and x,y G M with x 7^ y, we 
choose a minimal unit-speed g(t)-geodesic ^xy '■ [0,d g ^(x,y)] — > M from x to y. Note 
that we can choose 7^ so that (x, y) 1— >■ 7^ is measurable in an appropriate sense (see 
e.g. [25]). We use the same symbol 7^ for its range r fxy{^,d g {t){x,y)\). 
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3.1 Uniform bound for escape probability 

The goal of this subsection is to show the following: 

Proposition 3.4 Hindoo lim sup Q ^ P [<tr < T 2 ] = 0. 

For the proof, we will establish a discrete analogue of a comparison argument for the 
radial process as discussed in [18]. In this subsection, we fix R > 1 sufficiently large so 
that dg^T^io, Xq) < R — 1 until the final line of the proof of Proposition 3.4. We also 
fix a relatively compact open set Mq C M satisfying (2.1). Set tq := fo A (1/2), where 
^o — fo{Mo) is as in Lemma 2.6. 

The first step for proving Proposition 3.4 is to show a difference inequality for the 
radial process d g ( t )(o,X a (t)) (Lemma 3.7). It will play a role of the Ito formula for the 
radial process in our discrete setting. We introduce some notations to discuss how to 
avoid the singularity of d g (.)(o, •) on {o} U Cut s (.)(o). For r > 0, let us define a set A' r , A" 
and A r as follows: 

A'; := {(t,x,y) G [T ± ,T 2 ] xM„xM„| d g{t) (x,y) > r} , 

J\ T := J\ T I 1 J\ r . 

Note that A r is compact and that d g r.)(-, •) is smooth on A r . For t G [Ti,?^] an d p G M, 
let us define o p G M by 



'Tg(y) iHt,o, P )tA' n 

otherwise. 



Q (t) . = ) 'OP \0.) \*> U >*'J * ~ro> 



For simplicity of notations, we denote o n (a) by o n . Similarly, we use the symbol 7 n for 

7 n (a) throughout this section. Note that (t,o p ,p) <£ Cut ST holds. Furthermore, it is 
uniformly separated from CutsT in the following sense: 

Lemma 3.5 There exist r\ > and S\ > such that the following holds: Let to,t G 
[Ti, T2] with t — to G [0, #1]. Let po G -B#_i(o) owrf p E B 5 ° (po)- Then we have 

(i) d g (t)(o,p) < e K( * _ * o) (d ff (to)(o,po) + d fl (to)(po,p)), 

(ii) (t,og, o) ,p) G A ri w/ienp £ fiS o) (o). 

Here k = k(Mq) > is given according to Lemma 2.2. 

By applying Lemma 3.5 to X a , we obtain the following: 

Corollary 3.6 There exist ao > and h : [0, «o] - ^ [0, 1] u^^/i lim^o h(ot) = snc/i t/iat 
the following holds: For a < a , n G N and s,t G [t„ , t„ Q ji 1 ] 7 w/ien n < ctr, 



;i) d 9(t) (o,X a (s)) < e-* 2 (rf 9(t ( a)) (o,X°(ti a) )) + fc(a)), 
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(ii) (t,o n ,X a (s)) E A ri when X a {t { * ] ) £ B^ )] (o). 
Here 77 is the same as in Lemma 3.5. 

Proof of Corollary 3.6. Set Z := s\xp te [Ti,T 2 ],xeM 1^(0 L(t) ( rc )- Note that we have 
d t A a) AX a (t { n ) ),X a (t)) < y/m + 2a + Za 2 by the definition of X a . Take a > so that 
\/m + 2ao + Zaf t < Si and a 2 < 5± hold, where Si is as in Lemma 3.5. Then the conclusion 
follows by applying Lemma 3.5 with t = tn , Po = X a (tn ) and p = X a (s). □ 

Proof of Lemma 3.5. We show that (i) holds with S\ = 1. By the triangle inequality, 
the proof is reduced to showing the following two inequalities: 

d g{ t)(o, Po ) <e^- to U g{to) (o,p ), (3.1) 

d g( t)(Po,p)<e K ^ to) d g(to) (p ,p). (3.2) 

Our condition (2.1) yields that 7o P ° is included in Mo- Thus Lemma 2.2 yields (3.1). 
When p E B{ { (po), we have p E B R ° (o). Hence (2.1) and Lemma 2.2 yield (3.2) in a 
similar way as (3.1). 

Let us turn to consider (ii). For simplicity of notations, we denote o P q by d in this 
proof. We assume that t — t E [0,S] and p E Bg'(p ) hold for S > 0. First we will 
show (£, o',p) E A" , 4 when 5 is sufficiently small. Note that (t , o',p ) E A" , 2 holds since 

Po <£ -Bro (o) and d g ( to )(o,d) E {r /2,0}. Let q E 7^. By the triangle inequality, 

d g ( t )(o,q) < d g (t)(o,o') +d g (t)(o',po). (3.3) 

Since r /2 < 1 < i? holds, (2.1) yields 7^, C M when o' 7^ o. We can easily see that 
%,° po C 7op° C M . Thus, by applying Lemma 2.2 to (3.3), 

d B (t)(o,q) < e K(t ' to) (d g{to) (o,d) + dg {to) (d,p )) 

= {R-l)e Ka \ (3.4) 

Take 5 2 := 1 A (/sT 1 log(R/(R - 1))). Then, for any <5 G (0,5 2 ), (3.4) and (2.1) imply 
%, C Mq. Hence the triangle inequality, Lemma 2.2 and (3.2) yield 

dg(fi)(0',p) > dg(t)(0',P Q ) -dg(t)(pQ,p) 

> e-^- to U g{to) (d : p ) - e K ^d g(t0) (p ,p) 

> ^^ - e K5 S (3.5) 

when S < S 2 . Thus there exists S 3 = S 3 (n,r ,R) E (0,S 2 ] such that the right hand side 
of (3.5) is greater than r /4 whenever S E (O,^). Hence (t,o',p) E A" Q , 4 holds in such a 
case. 

Next we will show that there exists r[ > such that (t, o',p) E A' , holds for sufficiently 
small S. Once we have shown it, the conclusion holds with 77 = r[ A (r /4). As we did 



H:= \(t,x,y)e [f\,T 2 ] x M x M 
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in showing (t,o',p) G A" , 4 , we begin with studying the corresponding statement for 
(t , o',po). More precisely, we claim that there exists r'{ G (0, 1) such that (t , o',po) G A r » 
for all 5 G (0,1). When d = o, (to,o',po) G A' ro directly follows from the definition of 
d = off. When d ^ o, set 

r < d g{t) (o,y) <R-1, d 9 (t)(o,x) = r /2, 
d g (t)(x,y) = d g (jt)(o,y) - d g (t)(o,x) 

Note that H is compact and that H D CutsT — holds since (*, x,y) G H implies that 
x is on a minimal g(t)-geodesic from y to o. Since (*o,o',Po) G H by the definition of 
o', it suffices to show that there exists f\ > such that if C A~ . Indeed, the claim 
will be shown with r'[ — f\ l\r§ once we have proved it. Suppose that H C A' r does not 
hold for any r G (0,1). Then there are sequences (tj,Xj,yj) G if, (t'j , x'j , y'j) G CutsT, 
7 G N such that \tj — tj\ + d g (t j )(xj,x'j) + d g (tj)(yj,yj) — > as j — > oo. We may assume 
that ((tj,Xj,yj))j converges. Since (tj,Xj,yj) G H, x'j,y'j G M holds for sufficiently large 
j. Thus we can take a convergent subsequence of ((£'-, x'-,y'-))_j. Since CutsT and if are 
closed and d g (.)(-, •) is continuous, it contradicts with if n CutsT = 0- 

To complete the proof, we show that there exists Si G (0, £3] such that (t,o',p) G 
-4J.///2 when 5 G (0,<5i). Suppose that there exists (t',x',y') G CutsT such that |* — t'\ + 



d g (t)(o', x') + d g (t){p, y') < r i/2- For any q G 7^ ,, the triangle inequality and the assertion 



d g(t) (o,q) <d g{t) (o,p) + d g(t) (p,y') < e K5 (R - 1 + 5) + r'(/2. (3.6) 

A similar observation implies d g ( t )(o,q') < (e K< V + r'{)/2 for d G %, x ,- Thus there is 
64 = S^k, R) G (0, S3] such that the right hand side of (3.6) is less than R and (e K< V + 
fi)/2 < i? whenever 5 G (0,<5 4 ). In such a case, 7^ , C M and 7^, C M hold. Since 
(to,o',p ) G A^.„, Lemma 2.2 yields 

\t-t'\+d g{ t){o\x') + d g(t) (p,y') 

> |i - *'| - ^ + e- K5 rf 9(to) (o', a/) + e- K6 d g{to) (p, y') 

> e~ KS r'( + (1 - e~ KS )\t -t'\-5- e~ KS 5 (3.7) 

Take 5\ = Si(k, r") G (0, 5 4] so that the right hand side of (3.7) is greater than r'//2 when 
5 G (0,Si). Then (3.7) is absurd for any 5 G (0, Si). Thus it implies the conclusion. □ 

We prepare some notations for the second variational formula for the arclength. Let 
V^ be the g(£)-Levi-Civita connection and H® the g(t)-curvature tensor associated with 
y (*) _ p or a smooth curve 7 and smooth vector fields U, V along 7, the index form fy (U, V) 
is given by 

I?(U,V) := f ((V®U,V®V) m - (RV(U,^,V) g{t) ) ds. 

We write l\ (U, U) =: l\ (U) for simplicity of notations. Let G t , x ,y(u) be the solution to 
the following initial value problem on [0,d(x,y)\: 

r „ ( x _ ^C g(t) (jS(u),^( U )) 

G t , x , y (0) = 0, G' t>x>y (0) = 1. 



Coupling under a backward Ricci flow 13 



Note that G tjXjy (u) > for u G (0,d(x,y)] if y ^ Cut g ( t )(x) (see [18, Proof of Lemma 9]). 
For simplicity, we write G n := G («) V a^(«)v When X Q (t„ ) # _Br- ™ '(o), we define a 
vector field V^ along 7 n for each V G T , j a ),M by 

where /^™ 'V is the parallel vector field along 7 n of V associated with V^ tn \ Take 
f G M m . By using these notations, for n G No with n < N^ a \ let us define A n+ i and A n+ i 
by 

^n+l '■= (£n+l,7n) 9 (4"))' 

An+i := ^d ff(t ( i «») ) (o,o n ) + 9 t rf 9(t M ) (o n ,X Q (t^ ) )) 

+ (^(4 Q) ),7n) 9(t (.) ) (X Q (e))) + i/(f )(U +1 ) 

when X a (tn ) ^ -Br " (o), and A n+i = vm + 2(£ n+ i,?; )p and A n+1 = otherwise. 
Lemma 3.7 If n < ctrAN^', a < a® is sufficiently small and X a [tn ) ^ -Br ™ (o), t/ien 

rf 9(e+)i) (0,X°(tSl)) < ^(O^tt + <**«+! + « 2 A n+ l + 0(« 2 ) 

almost surely, where a^ is as in Corollary 3.6. In addition, o(a 2 ) is controlled uniformly. 

Proof. By virtue of Corollary 3.6, for sufficiently small a, the Taylor expansion together 
with the second variational formula yields 



d 



^WjK*"^)) < d g{t(na)) (o n ,X a (t^)) + a\ n+1 + a 2 d t d g{Aa)) (o n ,X a (t^)) 



+ a\Z{t^ n ) g{tia)) {X«{t^)) + pf\J in+1 ) 

+ o(a 2 ), (3.1 



O n ) = 



where J; 1 is a g(tn )-Jacobi field along 7 n with the boundary value condition J? 

and J? (X a (tn )) = £ n +i- Note that o(a 2 ) can be chosen uniformly since this expan- 
sion can be done on the compact set A ri and every geodesic variation is included in Mq. 

By the index lemma, we have J^J 1 (Jz J < !-/„ (Cn+i)- Hence the desired inequality 
follows when o n = o. In the case o n ^ o, we have 

d ff ( t w 1 )(°'- x ' a (*Si)) ^ ^(tl^)^' ") + rf 9 (4« + ) 1 )( «' XQ ( t i+i))' 

d^o, *"(#>)) = rf 9(t ( a)) (o,o n ) + d, ( #) ) (a n , *»(#>)). 

Note that (£« , o, o n ) is uniformly away from CutsT because of our choice of r and 
Lemma 2.6. Therefore the conclusion follows by combining them with (3.8). □ 
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Before turning into the next step, we show the following two complementary lemmas 
(Lemma 3.8 and Lemma 3.9) which provide a nice control of the second order term A n in 

Lemma 3.7. Set A n = E[A n | ^ n -i]- 

Lemma 3.8 Let (a n ) neNo be a uniformly bounded JP n -predictable process. Then 

N+l 



lim a 2 sup < \_, a j(Aj ~ Aj 



]=n 



n 



,N eN,n< N < N (a) A a R } = in probability. 



Proof. Note that the map (t,x,y) >->■ G t , x , y (d(x,y)) is continuous on A ri . Since we 
have G tjXj y(d(x,y)) > on A ri , there exists K > such that K~ x < G t;X;y (d(x,y)) < K. 
This fact together with Corollary 3.6 yields |Aj| and IAJ are uniformly bounded if j < 
ctr. Since X^=i a j(Aj ~~ A,-) * s an ^n-martingale and <7r is J^ ra -stopping time, the Doob 
inequality yields 



lim a 2 sup 



a-^-0 



N+l 
3=1 



in probability. 



(3.9) 



0<7V <N (°0 Aa R 

Here we used the fact linto^o^ 2 ^ — T 2 — 7\. Note that 



N( a ^Aa R N 

U U 

JV=1 n=l 



a 



N+l 



^djiAj-Aj 



]=n 



>5 



N^Aaa N 

c u u<- 

JV=1 n=2 



n-1 



a sup 

0<N<N( a )Aff R 



^ajiAj-Aj) 

=i 

N 

^o,(A,-A 



i=o 



> - >U <<r 



5 



N+l 



5 
>2 



Thus the conclusion follows from (3.9). 



□ 



Lemma 3.9 There exists a deterministic constant Co > being independent of a and R 
such that the following holds: 



♦(«)<> 



1 f d g(t (°<) ) ( > xa( - tn )) 
A n+ i < C + 77 / '"' " b{u)du. 

* Jo 

Proof. By using (m + 2)E[(£ n , ej)(£ n , e^)] = <5^j, we obtain 



E 



i=2 ^ 



(m-l)C;(d(o n ,X»(ti a) ))) 
G„(rf(o n ,X«(ti Q) ))) 
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Note that we have 

(z(tP),i n ). tM JX"(tP)) - (z(t^)n n ) n(f(a) Jo n ) 



' 9 {try 



' 9 (try 

d , (a) (o n ,X"{t ( n a) )) 



d^K.X"^)) _ u(a) 



^(^(4 a) )>7n) 9(i { i «) ) (7n(s))|s=«C?M 

(V(f ) Z(^)),7 n ) fl(t ( a)) ( 7n ( W ))^. 



Recall that, for (t,x,y) <£ CutsT, we have 

1 f d g(t)( x ,y) 

d t d g (t){x,y) 



2 Jo 



(dtg(t))(i%(u),^(u))du 



(cf. [19, Remark 6]). By combining them with Assumption 1, 



An+i = dtd^w^o, o n ) + - 



I r d u(a) (o n ,x«(t^)) 

J- / a(V„ ') 



%(4 Q) )(7n(M),7n(M))^ 






2G n (d(o n ,X*(C>))) 



A<*h 



<- / "' " ' b^du + dtd^^o, o n ) + (Z^),^)^)^) 



(«)> 



1 r d (a) ( O „,x«(4"0) 

1 / 9(*n ) 



+ 



2 Jo 
(m-l)G^(o n ,X*(d a) ))) 



Ric ff(t (a) ) (7„(«),7„(w))du 



2G„(rf(o n ,X-(d a) ))) 
Here we used the fact b(u) > in the case o n ^ o. Note that 

(m - l)G;(r) 



(3.10) 



Ric /.(a)x(7„(«),7„(«))d«- 



7/(C ; ) 



G„,r 



is non-increasing as a function of r. Indeed, we can easily verify it by taking a differenti- 
ation. Set 



/ 



C\ := sup sup 
tepri ) T 2 ] a . 6B (t) (o) 



\ 



Z(t)\ g(t) (x)+ sup (d t g(t)(V,V) + \mc g{t) (V,V)\) 



\ 



V£T X M 



/ 



By virtue of Lemma 2.2, C\ < oo holds. By applying a usual comparison argument to 
G' n (ro)/G n (r ), we obtain 



d (a) (o n ,X a (t ( n a) )) 

Sl*n ) 



Ric , j a )J'j n (u),'j n (u))du + 



{m-l)G> n {d{o n) X«{t^))) 



^C^ro + cotMdro)). 
Hence the conclusion follows from (3.10) with Co = C\(l + 3ro/4 + coth(Ciro)/2). □ 
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In the next step, we will introduce a comparison process to give a control of the radial 
process. Let us define two functions ip and ip on (2r , oo) by 

1 r 

if(r) :=C + - / b(u) du, 



ip{r) :- 



r — 2r, 



i 



where Cq is as in Lemma 3.9. Let us define a comparison process p a (t) taking values in 
[0, oo) inductively by 



p a (0) :=d 9 (o)(o,xo) + 3r , 

P a (t) := P a (t ] ) + 1 —^- («A n+1 + a\ V {p a {t^)) + ^(p a (t^)))) , t G [#>, *2il- 



a' 

The term ip(p a (tn )) is inserted to avoid a difficulty coming from the absence of the 
estimate in Lemma 3.7 on a neighborhood of o. By virtue of this extra term, p a (t) > 2r 
holds for all t G [Ti, T 2 ] if a is sufficiently small. Let &' R and a' R be given by a^, := cr R (p a ) 
and (T^, := CK~ 2 (|_5#Ja ~~ ^i) + 1- The following is a modification of an argument in the 
proof of [11, Theorem 3.5.3] into our discrete setting. 

Lemma 3.10 For 5 > 0, there exist a family of events (E£) a with \im a ^ Q ¥[Eg] = 1 and 
a constant K(S) > with lim^o K(S) = such that, on Ef, 

d g(t) (o,X a (t))<p a (t) + K(5) 

for t G [Ti, <7r A a' R A T 2 ] and sufficiently small a relative to 5 and R~ l . 

Proof. It suffices to show the assertion in the case t = £« for some n E Nq. Indeed, 
once we have shown it, Corollary 3.6 (i) yields 

d g{t) (o : X<*(t)) < e^ 2 {d gma) (o,X a ([t\ a )) + h(a)) 

< p « tla+ K(S) + (e Ka2 -l)R + e™ 2 h(a) 
< p a t + K(6) + a + (e Ka2 - l)R + e Ka2 h(a) 

for t G [Ti, &r A T 2 ]. Here we used the fact ip > and ^ > 0. Thus the conclusion can be 
easily deduced. 

For simplicity of notations, we denote d . .(<*), (o, X a (tn )) and p a (tn ') by <i n and p n 
respectively in the rest of this proof. Let us define a sequence of ^-stopping times Si by 
Sq := and 



S 2 i+i := inf <^ j > S 2 i 
S 2 i :=inf <j > S 2 i-i 



X a {& ) )eBif ) \o)\AN^ 



(«K 



X a (tf)t 3^,1(0)} A KM. 
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Since p n > 2r , it suffices to show the assertion in the case S 2 i < n < S 2 i+i A cr R A <j' r for 
some / G No- Now Lemma 3.7 and Lemma 3.9 imply 

d j+i ~ Pj+i < d j ~ Pj + <* 2 {}P(Aj) ~ <P(Pj)) + « 2 ( A i+i - Aj+i) + °(« 2 ) 

for j G [5*2/, S21+1 Ac^A er^). Here we used the fact ip > 0. Let f a be a C 2 -function on R 
satisfying 

VV ,/a|(— oo,— a) — U, Ja| (ckjOo)^/ •kj 

(ii) / Q is convex, 
(iii) a 2 sup x6R f£(x) = 0(1). 
For example, a function f a satisfying these conditions is constructed by setting 



fa) 



x pt 



bexp 



—00 j —00 



1-s 5 



l(-i,i)(s)dsdt, 



where a, b is chosen to satisfy 



cxp 



l-s< 



l(_i ; i)(s)ds = 1, 6 



1 r t 



cxp 



00 •/ —00 



1-s 5 



l(-l,l)(s)rfsrf£ = 1 



and f a (x) := af(a l x). By the Taylor expansion with the condition (iii) of f a , we have 

fa{d j+ i - pj+l) < f a (dj - Pj) 

+ a 2 f a (d 3 - P j){ip{dj) - ip{Pj) + (A, - A,)) + o(a 2 ). (3.11) 

Let C > be the Lipschitz constant of (p on [0, R]. Note that we have 

fM - Pj)(<p(dj) - <p(Pj)) < C(d 3 - Pj ) + + o(l). (3.12) 

Here the error term o(l) may appear in the case dj — pj G [— a,0]. Now by using (3.11) 
and (3.12) combined with the fact ds 2l — ps 2i < ~ a f° r sufficiently small a, we obtain 

\d n — Pn)+ < Ja\d n — Pn) 

n—1 n—l 

< Co 2 J2 (dj ~ Pjh + o? Y. f'M - ft)(Vi - A i+ i) + o(l). (3.13) 

j = S'2k j=S2k 

Here the first inequality follows from the condition (ii) of f a and n < a~ 2 (T 2 — T{) is used 
to derive the error term o(l). Let Ef be an event defined by 



Ex := < of sup 

k<k'<N( a )Aa R 



J^fM-i-Pj-iXAj-Aj 

j=k 



<s . 
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Note that a,j = f' a {dj-\ — Pj-i) is ^-predictable and uniformly bounded by 1. Thus, by 
combining Lemma 3.8 with (3.13), we obtain 

n— 1 
(d n - Pn)+ < CO? ^ ( d J ~ Pj)+ + 25 

j=s 2l 

on Ef for sufficiently small a. Thus, by virtue of a discrete Gronwall inequality (see [28] 
for instance), 

(dn - Pn)+ < 25 (1 + (1 + Ca 2 ) n ) < 25(1 + e c(T2 " Tl) ). 

This estimate implies the conclusion. □ 

Corollary 3.11 For every R' < R, lim sup a ^ F[a R < T 2 ] < lim sup Q ^ F[a' R , < T 2 }. 

Now we turn to the proof of our destination in this section. 

Proof of Proposition 3.4. By Corollary 3.11, the proof of Proposition 3.4 is reduced to 
estimate P[<t^. < T 2 ]. To obtain a useful bound of it, we would like to apply the invariance 
principle for p a . However, there is a technical difficulty coming from the unboundedness 
of the drift term of p a . To avoid it, we introduce an auxiliary process p a in the sequel. 

Let (p be a bounded, globally Lipschitz function on R such that <p(r) = (p(r) + ip(r) 
for r G [2ro + -R -1 , R]- Let us define an R- valued process p a (t) inductively by 



p a (0) := d g(Tl) (o,x ) + 3r , 

P a (t) ■= P a (t { n ] ) + 1 ~^f- («A n+1 + a 2 p(p a (t^))) , t e [t^^l] 



We also define two diffusion processes p°(t) and p°(r) as solutions to the following SDEs: 

dp°(t) = dB{t) + (<p(jP(t)) + iP(p°(t)))dt, 
p°(Ti) = d fl (Ti)(o,x ) + 3r , 
dp°(t) = dB{t) + ip(p°(t))dt, 
p°(Ti) = d g ( Tl )(o,x ) + 3r , 

where (Bit))^^^] i s a standard 1-dimensional Brownian motion with B(T{) = 0. We 
claim that p a converges in law to p° as a — > 0. Indeed, we can easily show the tightness 
of {p a ) a >o by modifying an argument for the invariance principle for i.i.d. sequences since 
(p is bounded. Then the claim follows from the same argument as we used in the proof of 
Theorem 3.1 under Proposition 3.2, which is based on the Poisson subordination and the 
uniqueness of the martingale problem. 

Let us define n R : ^ -»■ [7\, T 2 ]u{oo} by n R (w) = inf {t G [T U T 2 ] \ w(t) < 2r + R- 1 }. 
Then we have 

P[<^ < T 2 ] < P [a R (p a ) A VR (p a ) <T 2 ]=F [a R (p a ) A m {p a ) < T 2 ] . 

Since {w \ cr R (w) A n R (w) < T 2 } is closed in ^i, the Portmanteau theorem implies 

lim sup P [a R (p a ) A n R (p a ) < T 2 ] < P [a R (p°) A r] R (p ) < T 2 ] = P [a R (p°) A ^(p°) < T 2 ] . 



-,0 : c 



Since p is a diffusion process on (2r , oo) which cannot reach the boundary by Assump- 
tion 1, the conclusion follows. □ 
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3.2 Tightness of geodesic random walks 

Recall that we have metrized the path space ^ by using rf g m). To deal with the tightness 
of (X a ) ae (o,i) in ff, we show the following lemma, which provides a tightness criterion 
compatible with the time- dependent metric d g ( t y 



Lemma 3.12 (X a ) Q , G (o,i) is tight if 



lim - lim sup sup P 

<5^0 a-xx, n£ N 



sup d g{s) (X a (t^),X a (s))>e,a R 

tk a) <s<(t ( n a) +S)AT 2 



'X 



(3.14) 



holds for every e > and R > 1 . 



Proof. By following a standard argument (cf. [4, Theorem 7.3 and Theorem 7.4]), we 
can easily show that (X a ) Q , e ( 01 ) is tight if, for every e > 0, 



lim - lim sup sup P 

5^0 6 o^O ts[Ti,T 2 ] 



sup d g{Tl) (X a (t),X a (s))>e 

t<s<(t+6)/\T 2 



0. 



Thus, by virtue of Proposition 3.4, (X a ) ae ( 0j i) is tight if 



lim - lim sup sup P 

5^0 5 o^o tS[Ti,T 2 ] 



sup d g{Tl) (X a (t),X a (s)) > e,fr R = oo 

t<s<(t+S)AT 2 



for every £ > and R > 1. Given i? > 1, take Mo and k; as in Lemma 2.4 and Lemma 2.2 
respectively. Then, for e < 1 and s,t G [Ti,T 2 ], 

{rf 9(s) (X Q ( S ),X Q (LtJ Q )) < £,<r« = oo} C {rf 9(Tl) (X a ( S ),X a (t)) < 2e K(T2 - Tl) £,a« = oo} 

if a is sufficiently small. Thus we have 



sup d g(Tl) (X a (t),X a (s)) > e,a R = oo 

t<s<(t+<5)AT 2 



C 



sup 



[H a <s<([ti a +25)AT 2 

for a 2 < 5 and hence the conclusion follows. 



d 9 (s)(X a ([t\ a ),X a (s)) > G KT \ Tl£ ,a R = oo 



D 



Proof of Proposition 3.2. Take R > 1. By virtue of Lemma 3.12, it suffices to 
show (3.14). Take Mq C M compact and /t as in Lemma 2.4 and Lemma 2.2 respectively. 
By taking smaller e > 0, we may assume that e < ro/2, where fo = fo(Mo) is as in 
Lemma 2.6. Take n G No with n < N^ a ' . Let us define a ^-stopping time £ e by 



C £ := inf A; G N | n < k < N™, d ,, a) JX a (t^),X a (t^ } )) > 



9(trr 
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Then, for sufficiently small a, 



sup d g{s) (X a (t^),X a (s)) > 2e,a R = oo 

C {a 2 {( £ -n) <5, a R = oo}. (3.15) 



A a) <s<(tk a) +5)AT 2 



Set pk := X a (t k ) for fc G N and f(t,x) := d g ^(p n ,x). Note that / 2 is smooth on 
{/ < e}. Let us define \' k by 

\+l := \Sfc+l?7p».Pfc /gt^Y 

We claim that there exists a constant C > such that 

M+Wi) 2 < M^) 2 + 2«M a) , Pfc )A' fc+1 + ^« 2 (3-16) 

for k < Ce A JV( a ) on {<tr = oo}. Indeed, in the same way as we did to obtain (3.8), 

M+Wi) 2 < M°U) 2 + 2af(t k a \p k )X' k+1 + a 2 (X' k+1 ) 2 

+ 2a 2 f(t k a \p k ) (d t f(t k a \p k ) + (Z(ti a) ),^) g(ia)) (p k )) 

+ a 2 f(t k a \p k )I { fhj L+i ) + o(a 2 ). (3.17) 

'PnPk 

Here o(o; 2 ) is controlled uniformly. Let K\ > be a constant satisfying that the #(£)- 
sectional curvature on M is bounded below by —K 1 for every t G [7i, T 2 ]. Such a constant 
exists since Mo is compact. Then a comparison argument implies 

-y k 
'PnPk 

Here the right hand side is bounded uniformly if k < ( £ A N^ . The remaining estimate of 
the second order term in (3.17) to show (3.16) is easy since we are on the event {a R = oo}. 
Applying (3.16) repeatedly from k = n to k = ( £ , we obtain 

e 2 <aY,f(t k a \ Pk )K + i + CS 

k=n 

on {« 2 (Ce — n) < 5, &r = oo}. Set Ng := sup{£; G N | k < cr 2 b + n}. By taking 
5 < (2Cy 1 e 2 , we obtain 

{a 2 (( £ - n) < 5, a R = oo} C I ^2 f(t[ a) ,p k )\ k+1 > — , a 2 (( £ - n) < 5, a R = oo 



2a' 

. fc=n 



£ 2 



sup V/(t[ tt) )ft )l ,, Ai +1 > — \ (3.1 
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Set 



Y k , 



1 •= 



,(«) 



Vm + 2 



f(h 'P fc ) 1 |f^ a) .o,K e l^fe+l- 



{/(*r>fc)<er 



1 JV 



We can easily see that |Yj;| < 1 and ^2 k=n+1 Y k is J^/v- martingale. By [10, Theorem 1.6] 
with (3.18), we obtain 



[a 2 (C £ - n) < 5, & R = oo] < P 



7V+1 



sup 



S> 



fe > 



n<AT<AT (a) fc =n+ i 



2a\/m + 2 



< exp 



< exp 



AVm + 2(ae 2 + 2«Vm + 2(Nf } - n)) 



Wm + 2(ae 2 + 2^/m + 26) 
Hence (3.14) follows by combining this estimate with (3.15). 



□ 



4 Coupling by reflection 

For fcgl, let U at k be a 1-dimensional Ornstein-Uhlenbeck process defined as a solution 
to the following SDE: 



dU a , k (t) = --U atk (t)dt + 2dB(t), 
U a ,k(Ti) = a. 

More explicitly, U a ,k(t) = e~ k(t - Tl)/2 a + 2 f* e k ^-^/ 2 dB{s). Here B(t) is standard 1- 
dimensional Brownian motion as in the proof of Proposition 3.4. 



Theorem 4.1 Suppose 



2(VZ(t)) b + d t g(t) < Bic g , t) +kg(t) 



(4.1) 



holds for some k G R. Then, for each x±,X2 £ M 7 i/iere exists a coupling X(£) := 
(X L (t),X 2 (t)) of two ^-diffusion particles starting at (xi,^) satisfying 



P 



inf d fl(t) (X(t)) > 

Ti<t<T yw 



<P 



inf C/d-m^x! ,«!,),*(*) > 



X 



/or eac/i T G [Ti,T 2 ] ; where 



1 ^ 



XW := 



X 27T 



e fct -l 



-" 2 / 2 rf M , 



/3(t) := 



^g(Ti) (^1,^2) 

2 V //3(T-T 1 ) 



fc = 0. 
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Remark 4.2 (i) Given k G R ; a simple example satisfying (4.1) can be constructed 
from a solution g(t) to the Ricci flow d t g(t) = Ric^) by a scaling. That is, g{t) = 
e~ kt g(t) satisfies (4.1) (with equality) when Z(t) = 0. 

(ii) Our assumption (4.1) can be regarded as a natural extension of a lower Ricci curva- 
ture bound by k. Indeed, Bakry-Emery's curvature- dimension condition CD(k, oo) 
(see [2] for instance), which is a natural extension of a lower Ricci curvature bound, 
appears in (4.1) when both Z(t) and g(t) are independent oft. 

(iii) From the last item in this remark, when Z(t) = 0, one may expect that (4.1) works as 
an analogue of Bakry-Emery's CD (k,N) condition, which is equivalent to Ric 9 > k 
and dimM < iV when g(t) is independent oft, instead of CD (k, oo) since dimM = 
m < oo in our case. However, the following observation suggests us that we should 
be more careful: Let us consider (4.1) in the case k > and Z(t) = 0. When 
dtg(t) = 0, the Bonnet-Myers theorem tells us that the diameter of M is bounded 
and hence M is compact. When g(t) depends on t, it is no longer true. In fact, 
we can easily obtain a noncompact M satisfying (4.1) with k > by following an 
observation in the first item of this remark. On the other hand, the Bonnet-Myers 
theorem is known to hold under CD(k, N) when Z is of the form Vh in the time- 
homogeneous case (see [3, 22]). 

By following a standard argument, Theorem 4.1 implies the following estimate for a 
gradient of the diffusion semigroup: 

Corollary 4.3 Let ({X (t)} te [ TltT2 }, {F x } x€M ) be a Jf t - diffusion process with F x [X(Ti) = 
x] = 1. For any bounded measurable function f on M, let us define P t f by P t f(x) : = 
K x [f(X(t))}. Then, under the same assumption as in Theorem 4-1, we have 



lim sup 



P t f(x)-P t f(y) 



d g {Ti)(x,y) 



< /9 m m sup \f(*)-fV)l 

A/27rp(t — Ti) z,z'eM 

In particular, P t f is d g <T{)- globally Lipschitz continuous when f is bounded. 

Proof of Corollary 4.3. Let X = (X U X 2 ) be a coupling of ^-diffusions (X(t),F x ) 

and (X(t),F y ) given in Theorem 4.1. Let r* be the coupling time of X, i.e. r* := inf{£ G 
[Ti,T 2 ] | X(t) G D(M)}. Let us define a new coupling X* = (X*,X*) by 

x , ( (XW Ifr-X, 

1 (Xi(t),Xi(t)) otherwise. 

Since {r* > T} = {ini Tl < t < T d g{t) {X{t)) > 0}, Theorem 4.1 yields 

P t f(x) - PJ(y) = E[f(X*(t)) - f(X* 2 (t))} 

= e [(/(xr(t)) - f(x;(t))) i {T . >t} ] 

<P[r*>t] sup \f( z )-f( z ')\ 

z,z'eM 

( d g{Tl )(x,y) \ , 

\2y/P(t-Ti)J z,z'eM 
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Hence the assertion holds by dividing the both sides of the above inequality by d 9 (Ti)(x, y) 
and by letting y — y x after that. □ 

As we did in the last section, let ("f X y)x,yeM be a measurable family of unit-speed 
minimal g(t)-geodesics such that 7^ joins x and y. Without loss of generality, we may 
assume that 7^ is symmetric, that is, r yxy(d g (t)(x,y) — s) = 7^(5) holds. Let us define 
mi] '■ T y M^T y Mhj 

m®v :=v- 2{v,>y®) g{t) >y®(d g{ t)(x,y)). 

This is a reflection with respect to a hyperplane which is g(£)-perpendicular to jiy. Let 
us define m X y : T X M — y T y M by m X y '■= rrvxy o // {t) . Clearly m X y is a #(£)-isometry. As in 

Ixy 

the last section, let $**•' : M — y 0^{M) be a measurable section of the g(t)-orthonormal 
frame bundle ^ W (M) of M. Let us define two measurable maps <&f> : MxM ->■ ff^\M) 
for % — 1 , 2 by 

$f } (x,y) :=$W(x), 

(t) = (mg^foj/), (x,j/)eMxM\D(M), 

2 ^,yj. ^(t)^), (x,y)eD(M). 

Take 07, X2 G M. By using $' • , we define a coupled geodesic random walk X a (t) = 



C +1 := V^+Mp* (X a (#>)) Cn+l, 



(Xf (t),X2(t)) by Xf (0) = X, and, for i G [CV n+1J , 

t («) 



for i = 1,2. We can easily verify that Xf has the same law as X a with x = Xi. 

In what follows, we assume (4.1). We can easily verify that it implies Assumption 1. 
Thus, by Theorem 3.1, (X a ) a> o is tight under Assumption 1. In addition, a subsequential 
limit X Qfc — y X = (X\, X2) in law exists and it is a coupling of two Jzfr diffusion processes 
starting at X\ and X2 respectively. We fix such a subsequence (a^ken- in the rest of this 
paper, we use the same symbol X a for the subsequence X Qfe and the term u a — y 0" always 
means the subsequential limit "a*; —y 0". Set a R := crR(d g ^(o,X"(-))) for i — 1,2. We 
fix i? > 1 sufficiently large until the beginning of the proof of Theorem 4.1. Let M C M 
be a relatively compact open set satisfying (2.1) for 2R instead of R. 

We first show a difference inequality of rf 9 ( t )(X Q (t)). To describe it, we will introduce 

several notations as in the last section. For simplicity, let us denote 7 n ',„■, ,„■, by 7 n . 
Let us define a vector field V n +\ along % by 

K+i := it )} (g +1 - (a+i,7n) ff( 4«.) ) 7n(0)) . 
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Take v G M m . Let us define A* +1 and A* +1 by 
A 



Vi+1 •" 



2<S+i, %>,(#)) if(yi,y 2 )^(M), 



2VmT2(^ n+ i,t)) otherwise, 



A 



1 



n+l •" 



d (a) ,(X«(tl a) )) 



(8#(#>) + 2(VZ(e))) b ) (%(*),%(«)) cfe 



+ 4!" a)) (K+o 



Tn 



1 {X-(il Q) )^D(M)}' 



For 5 > 0, let us define r 5 : ^ -> {T U T 2 } U {oo} by r 5 (w) := inf {t > 7\ \ w(t) < 5}. We 
also define f s by f s := r s {d g ^(X a (•))). 

Lemma 4.4 For n G Nq wit/i n < A r ^ Q \ we /lave 



kt 



(a) 



e"W 2 d 



+ e^ Q) / 2 («A; +1 + « 2 A; +1 ) + o(a 2 ) 



(4.2) 

w/ien n < fj A a R A o" 2 ^ and a is sufficiently small. Moreover, we can control the error 
term o(a 2 ) uniformly in the position o/X Q . 

Proof. When (t„ , X a (t„ )) ^ Cutsx, (4.2) is just a consequence of the second variational 
formula for the distance function combined with the index lemma for Ii, n . To include 

In 

the case (tn , X Q (tn )) G CutsT and to obtain a uniform control of o(a), we extend this 
argument. Let us define H and pi,P2 '■ H — ¥ [Ti,T2J x Mq x Mq by 

# := < (*, x, y,z) d g (t)(x,y)>8, 

( d g (t) (x, y) = 2d 9 (t) (x, z) = 2d g{t) (y, z) 

Pi(t,x,y,z) := (t,x,z), 

p 2 (t,x,y,z) := (t,y,z). 

If q = (t,x,y,z) G H, then Pi(q),p 2 (q) <£ CutsT since 2 is on a midpoint of a minimal 
g(£)-geodesic joining x,y. Since H is compact, Pi(H) and P2{H) are also compact. Hence 
there is a constant rj > such that 

(t,x,y)Ep 1 (H)Up 2 (H), v > 



inf ^ |t - f | H- d g(t) (*, *') + d s(i) (y, </) ( ^ ^ ^ £ ^^ 
Take a > sufficiently small relative to 77 and 8. Set 



Pr, 



7« 



^(t^)^"^"* )) 



(a) 



Pn := expJ£ K+i 



I 2 
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By the triangle inequality, we have 

^( XQ (4 Q) )) = d^y, (Xnt ] ),Pn)+d (tia)) (p n ,* 2 Q (4 Q) )) , 






(a) - 

+1- 



Since (£« , X a (t„ ),7 n (d , (a)s(X a (£„ ))/2)) G //", we can apply the second variational 
formula to each term on the right hand side of the above inequality Hence we obtain 
(4.2). For a uniform control of the error term, we remark that 7„ is included in M and 
the g(tn )-length of 7„ is bigger than 5. These facts follows from n < f$ A a R A a R and 
the choice of Mq. Thus the every calculation of the second variational formula above is 
done on a compact subset of [Ti,?^] x Mq x Mq which is uniformly away from CutsT- It 
yields the desired result. □ 

Let us define a continuous stochastic process U£ on M. starting at a by 

(n (a) 

We next show the following comparison theorem for the distance process of coupled 
geodesic random walks. 

Lemma 4.5 For each e > 0, there exists a family of events (E£) a such that P [Ef] 
converges to 1 as a — >■ and 

d g{t) (X a (t))<U^ Ti) ^ a{Ti)) (t)+e (4.3) 

for all t E [Ti , T2 A f,5 A o"^ A <r^] on £"" /or sufficiently small a. 



Proof. In a similar way as in the proof of Lemma 3.10, we can complete the proof once 
we have found E" on which (4.3) holds when t = tn € [Ti,T 2 A ?s A a R A a R \. Set 

A; +1 := E[A* +1 1 & n ] and A^ := 0. Then YT j= i e fc *^- l/2 (A* - A*) is an #"„-local martingale. 
Indeed, A* +1 is bounded if n < fr R A a R and so is A* +1 . Let us define E" by 



N+l 

+ («) 



£?:=< sup ^e«;^(A--A-)<^ 



4 a) <^A^Aa| 



i=i 



In a similar way as in Lemma 3.8 or [16, Lemma 6], we can show lim a _>o P [Ef] = 1. Since 
we have (m + 2)E[(£j, e fc )(£j, e/)] = Ski, we obtain 

a; + i < -|^ ( ^)(x a (e ) )). 



(a) 



Thus an iteration of Lemma 4.4 implies (4.3) on E° when t = t n ■ □ 
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Proof of Theorem 4.1. Take e G (0, 1) arbitrary. Let R > 1 be sufficiently large so 
that 

lim sup P [Zr Aa 2 R < T 2 ] < e. 



a^O 



It is possible by Proposition 3.4. Set a := d g (r 1 )(xi,X2)- Take T G P~i,T 2 ] and let 5 > 
be 5 > 2e. Then Lemma 4.5 yields 

P [t s > T] < P [{f s > T} n £ £ a n {4 A^> T}] + 2e 
< P [n /2 {U:) > T] + 2e. 



Thus we obtain 



inf U"(t) > 5/2 

te[T u T] - 



¥[t 5 >T]<¥ 

by letting e I 0. Note that L7" converges in law to U a as a — >■ 0. Since 
{w G CdTi,^] 4MxI) r«(d p( .)(w(.))) > T} 

is open and {w \ mi te [ Tlt T 2 ]w(t) > 5/2} is closed in C([0,T] — > R), 
theorem yields 



the Portmanteau 



P 



Ti in/ <T ^)(X(t))>5 



< liminfP[f 5 > T] 

a— >0 



< lim sup P 



inf U"(t) > 5/2 

te[Ti,T] 



< P 



inf UJt) > 5/2 

i€[Ti,T] 



Therefore the conclusion follows by letting 5 4-0. 



□ 



We can also construct a coupling by parallel transport by following our manner. In 
the construction of the coupling by reflection, we used a map m X y ■ By following the same 
argument after replacing m X y with // (t) , we obtain a coupling by parallel transport. The 

difference of it from the coupling by reflection is the absence of the term corresponding 
to A*, which comes from the first variation of arclength. As a result, we can show the 
following (cf. [16]): 

Theorem 4.6 Assume (4.1). For Xi,x 2 G M, there is a coupling X(t) = (Xi(t),X2(t)) 
of two J^f-diffusion particles starting at x\ and x 2 at time Ti respectively such that 

d g{t) (Mt))<e-^-^d g(s) (Ms)) 

for T\ < s <t < T 2 almost surely. 

It recovers a part of results studied in [1] . In particular, a contraction type estimate for 
Wasserstein distances under the heat flow follows. 



Proof. Let us construct a coupling by parallel transport of geodesic random walks 
X a = (X^Xg) starting at (xi,x 2 ) G M x M by following the procedure stated just 



Coupling under a backward Ricci flow 27 

before Theorem 4.6. By taking a subsequence, we may assume that X a converges in law 
as a — > 0. We denote the limit by X = (X 1 ,X 2 ). In what follows, we prove 



P 



sup (e fct /% (t) (X(£)) - e fc ^ (s) (X( S ))) > e 

T 1 <s<t<T 2 







for any e > 0. By virtue of the Portmanteau theorem together with Proposition 3.4, it 
suffices to show 



limP 



sup (e kt / 2 d m (X a (t)) - e ks ^ 2 d g{s) (X a (s))) >e,a 1 R Aoj t = oo 

Ti<s<t<T 2 



(4.4) 



for any R > 1. For simplicity of notations, we write d n := e ktn l 2 d .(<*).. (X a (£„ ')) in this 
proof. For 5 > 0, let us define a sequence of ^-stopping times Si by S :— and 

S 2l+1 := ini {j > S 2l \dj<S}AN^\ 
S21 := inf {j > S21-1 I d,- > 25} A JVK 



Note that ds 2l _ x < 35 holds on {o^, A o - ^ = 00} for sufficiently small a. As mentioned just 
before Theorem 4.6, Lemma 4.4 holds with A* = 0. Moreover, we can obtain the same 
estimate (4.2) even when S 2 i-i < n < S 2 i A a R A a R for some / G N . In this case, the 
error term o(a 2 ) is controlled uniformly also in I. Let us define an event Ef by 

N+l - 

E?:={ SU p ^e^(A'-A-)<A 

nKNKN^ j= n +l 

Then, as in Lemma 3.8 and Lemma 4.5, we can show lim a _>o Pt-E'j] — 1- On Ef fl {a R A 
a R = 00}, we have d^ < d n + 5 for S21-1 < n < N < S 2 i if a is sufficiently small. 
Moreover, for n < S21-1 < N < S 2 i, 

d N - d n < (d N - ds 2l _ x ) + ds 2l _ L < 55 

In the case S 2 i < N < S 2 i+i, we obtain d^ — d n < 25. Thus d N — d n < 55 holds for all 
n < N on Ef fl {& R Aa R — 00}. Take 5 > less than e/10. Then our observations yield 
(4.4) since d 9 (t)(X a (£)) — d g (\t} a )(X-(\t\ a )) becomes uniformly small on {a R Aa R — 00} as 

a ->■ 0. " ~ D 
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